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ABSTRACT 

We report on the first X-ray proper- mot ion measurements of the nont her mally- dominated forward 
shock in the northeastern limb of SN 1006, based on two Chandra observations taken in 2000 and 
2008. We find that the proper motion of the forward shock is about 0.48" yr _1 and does not vary 
around the rim within the ~10% measurement uncertainties. The proper motion measured is con- 
sistent with that determined by the previous radio observations. The mean expansion index of the 
forward shock is calculated to be ~0.54 which matches the value expected based on an evolutionary 
model of a Type la supernova with either a power-law or an exponential ejecta density profile. As- 
suming pressure equilibrium around the periphery from the thermally-dominated northwestern rim to 
the nonthermally-dominated northeastern rim, we estimate the ambient density to the northeast of 
SN 1006 to be -0.085cm- 3 . 

Subject headings: ISM: individual (SN 1006) — shock waves — supernova remnants — X-rays: ISM 


1. INTRODUCTION 

SN 1006 is a Galactic shell-type supernova remnant 
(SNR) originating from a Type la SN explosion. In 
X-rays, its northeastern (NE) and southwestern (SW) 
limbs are dominated by nonthermal synchrotron emis- 
sion. These provided the first solid evidence that shells 
of SNRs are efficient accelerators of cosmic rays (Koyama 
et al. 1995; Winkler & Long 1997b). After the discov- 
ery of the nonthermal synchrotron emission from the NE 
and SW limbs of SN 1006, TeV 7 - ray emission had been 
expected to be detected, but has not yet been estab- 
lished 7 (Aharonian et al. 2005). TeV emission associ- 
ated with SNR shocks can arise from either of two pro- 
cesses: inverse Compton scattering of cosmic microwave 
background photons by relativistic electrons (the lep- 
tonic process), or 7r°-decay emission resulting from ac- 
celerated protons interacting with the ambient medium 
(the hadronic process) . The TeV flux from the hadronic 
process strongly depends on the ambient density (Kseno- 
fontov et al. 2005). Thus, measuring the ambient density 
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is important for understanding the origin of any TeV 7 - 
ray emission detected from SN 1006. 

The ambient density around the thermally emitting 
NW and SE rim can be measured based on spectral mod- 
elling of the thermally-dominated spectra. On the other 
hand, for the NE and SW rims where nonthermal emis- 
sion dominates, it is quite difficult to measure the ambi- 
ent density from spectral analyses. We may, however, in- 
fer the ambient density around nonthermally-dominated 
rims indirectly by assuming constancy of where 

no is the ambient density and v s is the shock veloc- 
ity (i.e., a uniform pressure around the periphery from 
the thermally- dominated NW rim to the nonthermally- 
dominated NE rim). In this context, measuring the 
proper motion of the forward shock is important for 
estimating the ambient density around nonthermally- 
dominated rims of SNRs and constraining the origin of 
TeV 7 -ray emission. 

So far, proper-motion measurements of the SN 1006 
have been performed by optical and radio observa- 
tions. Optical proper-motion of the bright NW Ha fil- 
ament was first measured by Long, Blair, & van den 
Bergh (1988) and later by Winkler, Gupta, & Long 
(2003), who precisely determined the proper motion to 
be 0.280" ±0. 008" yr -1 along the entire length of the NW 
filament. The radio observations measured the mean 
proper motion for the entire rim of the remnant to be 
0.44"±0.13"yr~ 1 (Moffett et al. 1993). Here, we report 
on the first X-ray proper-motion measurements of NE 
rim of SN 1006, based on two Chandra observations sep- 
arated by about 8 yr. 

2. OBSERVATIONS 

We observed the NE rim of SN 1006 twice: first on 
2000 July 10 (ObsID. 732; PI: Long) and a second time 
on 2008 June 24 (ObsID. 9107; PI: Petre). The second 
observation was specifically intended to allow a proper 
motion measurement; we requested the same pointing di- 
rection, roll angle and exposure time as the previous ob- 
servation. The time difference between the two observa- 
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Fig. 1.- — Chandra two-color image after vignetting effects are 
corrected. Red and blue represent 0.3-1.0keV and 1. 0-8.0 keV X- 
rays, respectively. The image is binned by 1.97" and has been 
smoothed by a Gaussian kernel of a = 5.90". The intensity scale 
is square root. We measure proper motions in the white and yel- 
low box regions shown in this figure. The four reference stars are 
indicated as arrows (from PI to P4). 

tions was 8.04 yr. We start our analysis from level 2 event 
files processed with calibration data files in CALDB ver. 
3.4.0 for ObsID. 732, ver. 3.4.5 for ObsID. 9107. For- 
tunately, we see no background flares in the data from 
ObsID. 9107 so that we reject no data from the level 2 
event file for this data set. On the other hand, we exclude 
high-background periods for data from ObsID. 732. The 
resulting net exposure times for ObsID. 732 and 9107 are 
55.3 ks and 68.7ks, respectively. 

3. ANALYSIS AND RESULTS 

Figure 1 shows a two-color image obtained from the 
2008 observation, after exposure correction to remove 
vignetting effects. Red and blue correspond to a low- 
energy (0.3-1.0keV) and a high-energy (1. 0-8.0 keV) 
band, respectively. The nonthermal forward shocks ap- 
pear as bluish filamentary structures along the edge of 
the remnant. In the following analysis, we measure 
proper motions of various portions of the forward shocks 
by comparing the images obtained in the two epochs. We 
use an energy range of 1.0-8.0keV to measure proper mo- 
tions of the forward shocks. 

Although the absolute astrometric accuracy for the 
Advanced CCD Imaging Spectrometer reported by the 
Chandra calibration team 8 reached a good level of 90% 
uncertainty = 0.6", we further align the images taken in 
different two epochs by aligning positions of four point 
sources which are visible in the FOV. These four sources 
are marked as P1-P4 in Fig. 1. We determine their posi- 
tions at each epoch, using wavdetect software included 
in CIAO ver. 4.0. Three of them (i.e., PI, P3, and 
P4 marked in Fig. 1) have optical counterparts, based 

8 http://cxc.harvard.edu/cal/ 



Fig. 2. — Difference image between the first- and second- epoch 
in 1.0-8.0keV band. The image is binned by 0.492" and has been 
smoothed using a Gaussian kernel of a = 1.467". The intensity 
is linearly scaled from -0.5 to +0.5 counts bin -1 . Example one- 
dimensional profiles from the box area are shown in Fig. 3. 

on the Naval Observatory Merged Astrometric Dataset 
(NOMAD) catalog (Zacharias et al. 2005). The proper 
motions of the three sources are negligible, at less than 
~0.005"yr -1 (Note that the P3 source is identified as a 
QSO; Winkler & Long 1997a). Although we cannot find 
an optical counterpart for the P2 source, the X-ray posi- 
tions determined for the two epochs are consistent within 
the errors. Thus, we assume that the proper motion of 
the P2 source is also negligible. Then, we align the coor- 
dinates of the images so that the four point sources have 
the same coordinates in the two epochs, by using ccmap 
and ccsetwcs in IRAF IMCOORDS. In the alignment 
procedure, four parameters for the second-epoch image 
are allowed to vary: x and y shifts, rotation, and pixel 
scale. The best-fit parameters show very small differ- 
ences; both x and y shifts are less than 0.2", rotation 
is less than 50", and no modification is required for the 
pixel scale. After the alignment, we find the rms residu- 
als for the positions of the four point sources between the 
two epochs to be ~0.14". We take this value as a sys- 
tematic uncertainty in our proper- motion measurements. 

Figure 2 shows the difference image between the two 
epochs in the 1.0-8.0keV band, after being exposure- 
corrected, registered and normalized to match the count 
rates in the two epochs. The narrow black lines along 
the periphery in this difference image clearly indicate the 
motion of the forward shock. In order to measure the 
proper motions accurately, we generate one-dimensional 
X-ray profiles along the apparent direction of motion — 
extracted from 21 areas (white and yellow boxes in Fig. 1) 
along the forward shock. Yellow boxes indicate inner 
shocks where there are essentially two shocks at the same 
azimuth angle (presumably a projection effect). In each 
case we have extracted the profiles in a direction per- 
pendicular to the forward shock, as determined from the 
local tangent to the rim taken at each of the extraction 
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Fig. 3. — Example one-dimensional profiles extracted from the 
area shown in Fig. 2. Data points with triangles and circles are 
responsible for 2008 and 2000 observations, respectively. 


boxes. Varying this tangent direction by several degrees 
does not significantly affect the results. 

To calculate the expansion rates and expansion in- 
dices, we estimated the center for the SNR shell by fit- 
ting a circle to the outer surface-brightness contour of 
the whole remnant based on a ROSAT RRl image: RA 
= 15 h 02 m 54 s .9 (J 2000 ), DEC = -41°56'08.9 // (J 2000 ), 
shown in Fig. 1 and Fig. 2. This center lies within ~ 40" 
of centers for the entire SN 1006 shell estimated from 
radio (Roger et al. 1988), X-ray (estimated by eye on a 
ROSAT PSPC image: Willingale et al. (1996), and opti- 
cal (Winkler, Gupta, & Long 2003) images 9 . The num- 
bers of photons from the two epochs are equalized in 
each area by scaling. Each one-dimensional X-ray profile 
is binned by 0.492". An example one-dimensional profile 
is shown in Fig. 3, in which we can clearly see the mo- 
tion of the shock. Then, shifting the first-epoch profile 
and calculating the % 2 -values from the difference between 
the two profiles at each shift position, we search for the 
best-matched shift position for each area. The reduced 
X 2 -values, which are nearly equal to x 2_ value divided by 
the number of bins in each area (here 69), range from 0.5 
to 1.2 at the best-matched shift positions. In this way, 
we obtain a proper motion for each area. A detailed ex- 
planation of our method of proper-motion measurements 
can be found in Katsuda et al. (2008). 

Figure 4 left shows both proper motions and velocities 
(for a distance of 2.2 kpc; Winkler, Gupta, & Long 2003) 
for all the areas as a function of azimuthal angle, and 
Fig. 4 right shows both expansion rates and expansion 
indices, m; R oc t m and m is derived from (the expansion 
rate) x (the age of the remnant, 998 yr for observations 
used here). The data points with black and red are ob- 
tained from the white and yellow boxes in Fig. 1, respec- 
tively. The errors quoted are statistical 1 <t uncertain- 
ties. We find that the proper motion around the forward 
shocks is almost invariant. Also notable is that there is 
no significant velocity difference between the inner shock 
and the outer shock at the azimuth angle around 50°- 
70°. The weighted mean value of the proper motions of 
the forward shocks is calculated to be ~0.48"yr -1 with 
the rms deviation of ~0.04"yr -1 . The expansion rates 


9 Winkler, Gupta, & Long (2003) also give a center of curvature 

for bright NW optical filaments where they measured the proper 

motions; this position is offset ~ 5' SE from the shell center. 


TABLE 1 

Expansion Index Measurements for Young SNRs 


SNR name 

Radio 

X-rays 

Cas A 

0 . 35 1 

0 . 73 1 ’ 2 ’ 3 

Tycho 

0 . 47 4 

0.71 s 

Kepler 

0 . 5 6 

0.93 7 -*• 0.35-0.8 8,9 

SN 1006 

0.48 10 

0.54 11 


Note. — 1: Koralesky et al. (1998), 2: Vink et al. (1998), 3: De- 
Laney & Ludnick (2003), 4: Reynoso et al. (1997), 5: Hughes (2000), 
6: Dickel et al. (1988), 7: Hughes (1999), 8: Katsuda et al. (2008), 9: 
Vink (2008), 10: Moffett et al. (1993), 11: This work 


and the expansion indices are also nearly constant at 
about 0.054 %yr -1 and 0.54, respectively. The statisti- 
cal error is about 10% (~0.05"yr -1 ) for each data point, 
so that it is about three times greater than the system- 
atic uncertainty of 0.017" yr -1 (= 0.14"/8.04 yr). Thus, 
the current expansion index of the NE rim is consider- 
ably less than if it were in a free expansion (m = 1 ), 
but significantly higher than the value for an interstellar 
medium (ISM) -dominated adiabatically expanding SNR 
(m = 0.4). 


4. DISCUSSION 

This paper presents the first X-ray proper-motion 
measurement of the forward shock along the NE 
rim of SN 1006. The measured proper motion of 
~0.48"±0.04"yr -1 does not vary significantly with az- 
imuth. This value is consistent with the mean expansion 
rate for the entire SN 1006 shell measured by radio ob- 
servations, 0.44"±0.13"yr -1 (Moffett et al. 1993). For 
other young SNRs, there have been problems reconcil- 
ing expansion index measurements from X-ray and radio 
wavelengths, with the X-ray rate seeming to be larger 
in Cas A, Kepler’s SNR and Tycho’s SNR as noted by 
e.g., Hwang et al. ( 2002 ) or DeLaney et al. (2004). Ta- 
ble 1 summarizes the radio and X-ray expansion indices 
for these SNRs. We note that recent Chandrarb&sed X- 
ray proper-motion measurements have reduced the ex- 
pansion index for Kepler, more in line with the radio 
value. Cas A is the only SNR whose X-ray expansion 
index as measured with Chandra clearly differs from its 
radio value. Tycho’s SNR has yet to be remeasured with 
Chandra , but should be as soon as possible. 

We here estimate the ambient density around the NE 
rim of SN 1006 and discuss its implication for TeV 7 - ray 
emission. Along the NW rim, the proper motion of the 
forward shock is well determined to be 0.280" ±0.008" 
yr -1 by optical observations (Long, Blair, & van den 
Bergh 1988; Winkler, Gupta, & Long 2003). The ratio of 
the forward shock velocities between the NE rim and the 
NW rim is estimated to be ~1.7. The ambient density 
ahead of the NW rim is measured to be 0.15-0.25 cm -3 
by X-ray spectroscopy (Long et al. 2003; Acero et al. 
2007) or 0.25-0.4 cm -3 by optical spectroscopy (Ray- 
mond et al. 2007). We here take 0.25^0 . 10 cm -3 as a con- 
servative ambient density value. Then, assuming pres- 
sure equilibrium around the periphery from the NW rim 
to the NE rim (i.e., constant no^ 2 ), we estimate the ambi- 
ent density ahead of the NE rim to be 0.085to!o35 cm -3 . 
This is consistent with the upper limit, 0 . 1 cm -3 , con- 
strained by an upper limit of TeV 7 -ray flux obtained 
by HESS (Aharonian et al. 2005) combined with a 7 -ray 
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Fig. 4. — Left : Proper motion (left axis) and velocitiy for a distance of 2.2 kpc (right axis) as a function of azimuthal angle (counterclockwise 
from the north). The black and red data points are obtained from the white and yellow boxes in Fig. 1, respectively. The weighted mean 
proper motion (with rms) is shown in the panel. Right'. Same as left but for expansion rate (left axis) and expansion index (right axis). 


emission model (Ksenofontov et al. 2005). 

We derived the mean expansion index of the NE for- 
ward shock to be ~0.54. Dwarkadas & Chevalier (1998) 
examined the evolution of Type la SNe for both power- 
law and exponential density profiles of ejecta. From their 
calculations, the evolution of the forward shock is scaled 
by the time, t f = 248(no/l cm -3 ) -1 / 3 yr. The age of 
SN1006, 998 yr, corresponds to 1.48-2.08 depending 
on a range of the ambient density of 0.05-0. 14 cm -3 . 
Then, the expansion index is expected to be 0.57-0.53 
(power-law) or 0.54-0.51 (exponential). Thus, we find 
that the measured value is in good agreement with that 
expected in the Type la SNe model with either a power- 
law or an exponential ejecta density profile. 

We note that there are a number of thermally- 
dominated knotty features in the remnant (see, Fig.l). 
The origin of these features are suggested to be either 
fingers formed by Rayleigh- Taylor instabilities acting at 
the contact discontinuity between the swept-up ISM and 
the ejecta (Long et al. 2003; Cassam-Chenai’ et al. 2007) 
or ejecta bullets (Long et al. 2003; Vink et al. 2003). It 
is also interesting to measure proper motions for these 
features, which is beyond the scope of this paper and 
remains as our future work. 


5. CONCLUSION 

We, for the first time, performed X-ray proper-motion 
measurements of the forward shock in the NE limb of 
SN 1006. The proper motion of the forward shock is 
derived to be ~0.48"yr _1 . This is consistent with the 
results from radio observations. The expansion index of 
the forward shock is calculated to be ~0.54, which sup- 
ports that SN 1006 is a remnant of a Type la supernova 
with either a power- law or an exponential ejecta density 
profile. We estimate the ambient density to the NE of 
SN 1006 to be ~0.085cm -3 . Our derived ambient den- 
sity of 0.085 cm -3 is below the limit of 0.1 cm -3 obtained 
by Ksenofontov et al. (2005) in modeling hadronic and 
leptonic TeV 7-ray emission subject to the HESS upper 
limit (Aharonian et al. 2005). 
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